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The application of the Schwinger-boson transformation to quantum Heisenberg magnets is briefly 
reviewed, beginning with the derivation of a rotationally invariant mean-field theory. The inclusion 
of Gaussian fluctuations is discussed in some detail for a general (non-Bravais) lattice, extending 
available results for simple lattices. Numerical results are presented for the ground-state energy of 
the Shastry-Sutherland model, and possible future refinements of the method are outlined. 
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I. INTRODUCTION 

The Schwinger-boson method is a versatile means of 
studying quantum spin systems || [| Q. It relies on 
the correspondence between spin and boson degrees of 
freedom defined by the Schwinger-boson transformation. 
The introduction of fictitious boson degrees of freedom 
permits the derivation of simple mean-field theories, with 
rotationally invariant order parameters which cannot be 
defined in terms of the original spins. Because the pro- 
cedure explicitly preserves rotational invariance, sponta- 
neous symmetry breaking need not be assumed, in con- 
trast to standard spin- wave approaches. The method is 
therefore capable of describing collinear, spiral, and spin- 
liquid phases within a unified framework, where broken- 
symmetry phases are represented by Bose-condensed 
phases of Schwinger bosons 0. At the mean-field level, 
early Schwinger-boson calculations produced quantita- 
tively disappointing results. However, small changes in 
the formalism greatly improve its accuracy without un- 
duly complicating its use ||. In this paper we briefly 
review the application of the Schwinger-boson transfor- 
mation to quantum Heisenberg magnets, beginning with 
the derivation of a rotationally invariant mean-field the- 
ory. We then show in some detail how Gaussian fluc- 
tuations around the mean-field solution can be included 
for arbitrary (non-Bravais) lattices. This generalizes the 
available results, which treat simple lattices. As an appli- 
cation, we examine the Shastry-Sutherland model, which 
has been the subject of renewed interest due to its exper- 
imental realization in SrCu2 (603)2 @- 



II. SCHWINGER-BOSON METHOD 

The Schwinger-boson transformation is an exact map- 
ping from spin to boson degrees of freedom, given by 
S i— * |(aT b') m cr mn (a b) n ; here m,n — 1,2, and cr is 
the vector of Pauli matrices. The original spin Hilbcrt 
space corresponds not to the full boson Hilbert space 
but rather to the (2S + l)-dimensional subspace in which 
Sn = <va + b'b — 2S = 0. Our goal is to evalu- 



ate the partition function, Z((3) = Tr 



phys ' 



" , where 



7i({S;}) i— > H.({ai,bi}) is an arbitrary spin Hamiltonian 
transformed to a boson Hamiltonian, and the trace is re- 
stricted to the physical subspace. We can replace the 
restricted trace with a trace over the full boson Hilbert 
space using an appropriate projection operator for each 
spin: 



Vi = J exp (^-iXiSn^ . 



(1) 



Then Z = Tr (J],; Vi)e~ m , which can be rewritten as a 
functional integral in the usual way. In particular, we 
obtain 

_ I XtyqPfoPA S[d T X] exp {-A^g, fa, A]) 

JVX S[d T \] ' 1 > 

where the Euclidean action is 
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dr 
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(3) 

The Lagrange multipliers {A;} have been promoted to 
time-dependent variables; the inclusion of <5[<9 T A] = 
Y[ i S(d T Xi) in the functional measure ensures that Z 
remains the same. 

The action (^) is invariant under local gauge transfor- 
mations parametrized by {9i(r)}: X — > X—d T 8, <fi — > e %e <\>. 
This gauge symmetry is associated with the conservation 
of total boson number on each site, which holds simply 
because the magnitude of each spin is constant. The 
gauge is partially fixed in Eq. by the condition that 
d T X = 0; the remaining, unfixed, transformations are 
those where 9 is time-independent. The action will often 
have additional symmetries, such as global SU(2) invari- 
ance; these depend entirely on the symmetries of the spin 
Hamiltonian. 

We specialize to the Heisenberg Hamiltonian, H = 
X)<ij> JijSi ■ Sj. Under the Schwinger-boson transfor- 
mation Si ■ S 3 i * ■.Bl j B i j:-A] :j Aij, where A tj = \{aibj - 
bihj), Bij — \{a\cbj+b\bj), and :<D: is the normal-ordered 
form of O. The bilinear operators Aij and B^ are in- 
variant under simultaneous SU(2) rotations of spins i 
and j, so they may take on vacuum expectation values 
even without spontaneous rotational symmetry breaking. 
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This makes them convenient order parameters on which 
to base a mean- field theory. These operators are not, 
however, invariant under the U(l) gauge transformations 
described earlier (since they couple physical and unphys- 
ical states), so such theories retain some gauge degrees 
of freedom. These degrees of freedom become impor- 
tant when the order parameters are allowed to fluctuate 
around their average values. 

We write the Hamiltonian in a form suitable for lattices 
with several spins per unit cell: 



X, u 7 



(4) 



Here x labels the unit cell, \i and v = 1,2,... A s it C s are 
site labels within a unit cell, and n — 1,2,... Abonds is 
a bond label. We then formulate a mean-field theory 
with order parameters {A xra ,i? xn } by introducing com- 
plex Hubbard-Stratonovich fields a and /3, linearly cou- 
pling them to A and B, and using them to integrate out 
the original boson fields <f>. This yields a new representa- 
tion of the partition function, 

= J VaVpVX 8[d T \]exp(--Aes[a, P, A]) { , 
(P) ~ JVaV(3VX (5[a T A]cxp(-^Hs[a, 0]) ' 

where the effective and Hubbard-Stratonovich actions are 
Aff = ^ HS +TrlnM+ / dr V -2iS^, (6a) 

f/3 

+p: n p xn ). (6b) 

The denominator of (|^) can be evaluated exactly, but it 
is convenient for normalization purposes to leave it as it 
stands. The dynamical matrix M[a, /3, A], defined in the 
Appendix, describes the propagation of Schwinger bosons 
coupled to the Hubbard-Stratonovich fields. 

As noted by Trumper et at, the effective action (|6a| ) 
has a gauge symmetry inherited from the boson action 
(||) @- r 4' ne gauge-fixing functional, currently given by 
<5[c? T A], can be chosen more or less arbitrarily without 
changing the value of Z\ this can be done with the 
Faddeev-Popov procedure. However, while the exact 
value of Z is gauge- invariant, the terms in a general per- 
turbative expansion of it are not. We will therefore retain 
the current gauge- fixing condition in what follows, but 
observe that this quantitatively affects our results, be- 
cause the saddle-point expansion we use is not controlled 
by any gauge-invariant small parameter. By contrast, 
mean-field theories which take as order parameters cither 
{Axn} or {B xn }, but not both, have natural SU(N) gen- 
eralizations in which the saddle-point expansion is con- 
trolled by 1/N g. 

Next we look for stationary points of ,4 ff with respect 
to its arguments. We consider only static and transla- 
tionally invariant solutions of the form a xn = a* n = a n , 



A™ = —PZn = Pn, an d A XAI = -ik p , where a, /?, and 
A are real, and we work in the zero-temperature limit. 
With these restrictions we can derive the desired rota- 
tionally invariant mean-field equations. Collectively they 
take the variational form 6E = 0, where the energy per 
unit cell is given by 



N 



-E 



(7) 



and the quasiparticle energies {^k^} are defined in the 
Appendix. This represents a set of 2Abonds + -^sites 
coupled nonlinear equations, which must be solved nu- 
merically for the mean-field parameters a^°\ and 
A^ ). (There are generally several gauge-equivalent sad- 
dle points, differing only in the signs of a and /?; which 
one of these is chosen as the mean-field solution is irrel- 
evant.) The resulting boson propagator C? (0) = (M^)- 1 
enters into the calculation of fluctuation corrections. The 
ground-state energy per unit cell can be rewritten as 

^ (0) =E n X((/3i 0) ) 2 -(al 0) ) 2 ). 



Finally, we include fluctuations around the saddle 
point in the usual way, by expanding the action in a 
power series: 



A cfi = 



V d 2 A cS 
2 d^d^ 



* + 0(* 3 



(8) 



(o) 



where the fluctuation fields ^ = (a*, a, /?*, (3) — 
(a*, a, (3*, (3)^ do not include A due to our choice of 
gauge. The matrix of second derivatives, which is block- 
diagonal in reciprocal space, determines the Gaussian 
correction to the ground-state energy: 

+ OG 



where 



D am , fc „(k,.) = J-Tr gW^™-gW. ^ 



K.oj,arn 



ku . bri 



(10) 

These expressions, together with the details provided in 
the Appendix, represent our main analytical results. We 
now discuss the numerical results for a particular lattice. 



III. SHASTRY-SUTHERLAND LATTICE 

The Shastry-Sutherland model is a frustrated 2D an- 
tiferromagnet with bonds of two different strengths, as 
shown in Fig. |l|. The relevant physical parameters are 
the ratio g = J2/J1 and the spin magnitude S. In the 
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FIG. 1: Shastry-Sutherland lattice. Single- and double-line 
bonds have strengths Ji and J2 respectively. The dashed 
square contains one unit cell. 



classical limit (S — > 00) the ground state has Neel or- 
der for g < 2 and helical long-range order otherwise. 
For finite S and sufficiently large g, the exact quantum 
ground state is known: it is fully dimerized, with the two 
spins on each J2-bond in the singlet state ||. (This is 
often referred to as a spin-liquid phase due to the lack 
of long-range order.) The energy of this state is always 
— 2S(S + 1) J% per unit cell. Moreover one expects the 
Neel order present at g — to persist up to some positive 
g c {S). We would like to clarify to what extent quantum 
fluctuations melt the classical spiral ordering. Our main 
interest is in the case S — 1/2. This is the physically 
relevant case, and the existence of a spiral phase is in the 
greatest doubt here. A previous Schwinger-boson calcu- 
lation at the mean-field level showed a spiral phase in the 
range 1.1 < g < 1.65 ylj. This result has not been con- 
firmed by exact diagonalization, since the spiral states 
are only energetically favorable on large lattices. High- 
order series expansions, starting from both g = 00 and 
g = 0, have been used to locate accurately the energetic 
crossover between the Neel and spin-liquid phases | pl| ; 
however, these calculations provide no evidence for or 
against an intermediate spiral phase. 

When describing helical phases it is common to use 
"twisted" periodic boundary conditions, where spins 
along one cluster boundary are rotated with respect to 
those on the opposite boundary |j, [l2|]. This procedure 
breaks rotational invariance and complicates the analysis 
somewhat, so we do not use it here. The price we pay is 
that helices with small pitch cannot be accommodated by 
small clusters with untwisted boundary conditions. We 



found that helical mean-field solutions exist for clusters 
with as few as 24 sites, suggesting that finite-size effects 
probably do not seriously alter our results on larger clus- 
ters. 

We used the Schwinger-boson method to calculate the 
ground state energy for S — 1/2. The mean-field equa- 
tions can be solved numerically for the full range of g. For 
small g, the solution with lowest energy is purely antifcr- 
romagnetic (a ^ 0, (3 = 0) on the strong Ji-bonds and 
purely ferromagnetic (/3 ^ 0, a = 0) on the weaker J2- 
bonds. This solution is unique and has ordering wavevec- 
tor Q = (tt, 7t); we identify it as the Neel state. When 
g is increased, new solutions with ordering wavevectors 
Q = (tt ± Sq, tt) and Q = (tt, tt ± Sq) appear through a 
bifurcation of the Neel state; we identify these solutions 
as helical states with pitch along the x-axis and y-axis 
respectively. The helical states are degenerate for square 
lattices. Finally, for sufficiently large g, the energy of the 
helical states exceeds that of the spin-liquid state, and 
the ground state becomes disordered. 

The results for a small cluster (16-site, 4x4) are shown 
in Fig. ||. The Neel state does not undergo any bifur- 
cation here, due to the small size of the cluster. The 
Schwinger-boson mean-field estimate exceeds the exact 
ground state energy by a small and roughly constant 
amount over the full range of g. A second-order per- 
turbative calculation around 5 = (that is, around the 
Neel state) yields similar results in the displayed range 
(although naturally perturbation theory is more accurate 
near g — 0). Finally, the inclusion of Gaussian fluctu- 
ations in the Schwinger-boson calculation improves the 
accuracy, but for g < 1.3 only: for larger g the results 
deteriorate rapidly. 

In Fig. [| we show the mean- field and Gaussian-order 
results for a larger, rectangular cluster (64-site, 16 x 4), 
for which exact results are unavailable (except for the 
spin-liquid state). In this case the Neel state bifurcates 
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FIG. 3: Ground-state energy, large cluster. As in the previous 
figure, the solid curves represent Schwinger-boson mean-field 
results, while the dashed curves are corrected for Gaussian 
fluctuations. These energies are given for both Neel and heli- 
cal phases, as indicated. The exact energy of the spin-liquid 
phase is also shown (dotted line). 



at g ~ 1.2 into a helical state with pitch along the x-axis 
(the lower dashed curve). The energy of the Neel state 
is corrected by fluctuations just as in the smaller cluster: 
it is shifted downward by a small amount until g ~ 1.3, 
after which the computed correction becomes large and 
positive. Extrapolating from the small-cluster results, it 
seems probable that the upswing is again spurious, and 
that the true energy is bracketed by the mean-field and 
Gaussian-order results. The Gaussian correction to the 
energy of the helical state is uniformly positive. 

We also checked the mean-field and Gaussian-order re- 
sults for still larger clusters, at selected values of g. The 
bifurcation point occurs at slightly smaller g on larger 
clusters, approaching g = 1.1 in the thermodynamic 
limit, as seen in Ref. pp[ . The nature of the correc- 
tions does not appear to change. Our conclusion is that 
the Shastry-Sutherland model with S — 1/2 has a helical 
ground state, over a range of g narrower than that sug- 
gested by mean-field Schwinger-boson theory. The range 
of stable wavevectors for spiral ordering is correspond- 
ingly reduced: we estimate that the helical state is stable 
for 5q < 7r/3, whereas the mean- field theory predicts sta- 
bility for Sq < 37r/8. These conclusions are necessarily 
tentative. 



gies. 

The results can be interpreted in two ways. It may 
be that the very good agreement between the mean-field 
and exact ground-state energies is simply fortuitous; the 
partition function may not be dominated by contribu- 
tions from the neighborhood of the mean-field solution. 
In that case no refinement of the present saddle-point ex- 
pansion will yield better results. This is a real possibility; 
however, since the Schwinger-boson mean-field approach 
also gives accurate results on a number of other systems, 
it seems unlikely. 

It also may be that we have failed to treat the gauge 
degrees of freedom properly. As noted earlier, the choice 
of gauge-fixing condition has a quantitative effect on the 
order-by-order results of the saddle-point expansion. |l3| 
We worked in a simple but otherwise arbitrary gauge. It 
would be preferable to control the expansion in a gauge- 
invariant way. This should be possible along the lines 
of Ref. B (and references therein) : allow the number of 
boson "flavors" to become large, say N, and generalize 
the fields A, B, and 5n in a natural way to include these 
new flavors. When this is done correctly, the effective 
action for a, /3, and A will be proportional to N; the 
saddle-point expansion will then yield an asymptotic se- 
ries in 1/N, each term of which is gauge-invariant. This 
program would put the present work on a sounder the- 
oretical footing, by determining in what limit its results 
become exact. The key step is to generalize A and B 
simultaneously; this remains a challenge for the future. 
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IV. REMARKS 

From a practical point of view, these results are not 
very satisfactory. Although the inclusion of Gaussian 
fluctuations improves the agreement between Schwinger- 
boson theory and exact diagonalization in some regions 
of parameter space, the improvement is not dramatic. 
Moreover, the correction is not uniformly good; it gives 
little reliable information about the region of interest, 
where the Neel and helical states have comparable ener- 



APPENDIX A: DETAILS OF CALCULATION 

This appendix contains details about some aspects of 
the calculation which, while extraneous to the main body 
of the paper, are included for the sake of completeness. 

The boson action (which was integrated out) has the 
form ((/>* fa) ■ M • (cj) a (j)l) T ; i.e., it is anomalous. The 
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dynamical matrix M[a, /?, A] has 2x2 submatrices 





-iwe 



M^'A^ + ^(F+) k(tl , k ^/ (-ilm/3„)j i 
-^(F~)k M ,kv (-Rc/SnO'z - iRea„(T y + ilma n a x ) , 



(Al) 

where all fields are Fourier-transformed, with arguments 
(k — k', uj — uj') suppressed, and 

The differing convergence factors e ±lew in the first term 
of M arise from the anomalous nature of the action. They 
must be taken into accountproperly to obtain the cor- 
rect mean-field energy, Eq. (0) , and to evaluate the Gaus- 
sian correction given by Eq. (^J). However, in the latter 
equation only the first (n = 1) term from the expansion 
Tr ln(l — D) = — J2 n — Tr D™ needs these factors to force 
convergence of the w-integral, and the contribution to 
£;(!) from this term can be shown to equal ^E^°\ For 
numerical work, then, it is simplest to rewrite the Gaus- 
sian correction as 



2 



duj 

47T 



(Tr D + Tr ln(l — D) ) ( A3 ) 

k 



The mean-field dynamical matrix is block-diagonal: 



-^M (0) (k,w) = iu}l®& z +A{k)®i + B(k)®a V7 (A4) 



where A and B are Hermitian matrices obtained from 
the general expression for M. This can be diagonalized 
using a nonunitary (generalized Bogoliubov) transforma- 
tion: there is a matrix S(k) — G(k)CS>l + i?(k)ig>(7 y such 
that <S T = (i <g> <T2)5 _1 (i (g) a z ) and 









-ioj + Lo kfi 



(A5) 

The quasiparticle energies {^k M } are the square roots of 
the eigenvalues of (A + B)(A — B); these eigenvalues will 
be real and positive when A is large enough, and the 
mean-field solutions must be sought in that region. The 
diagonalizing matrix and the quasiparticle energies are 
analytically tractable for one- and two-site unit cells; in 
any case they are easy to obtain numerically. The boson 
propagator is then 



<V( k ) 



(iuj + bJ-kv)' 







-iuj + Wkv) 



<5L(k). 



and disregard the convergence factors. 



(A6) 

Finally, the mean-field derivatives of M with respect 
to the fluctuation fields are field- and frequency- 
independent, and can be read directly from Eq. (|Al[). 
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